Atomic-Scale Dynamics of the Formation and Dissolution 
of Carbon Clusters in Si02 
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Oxidation of SiC produces S1O2 while CO is released. A 'reoxidation' step at lower temperatures 
is, however, necessary to produce high-quality Si02- This step is believed to cleanse the oxide of 
residual C without further oxidation of the SiC substrate. We report first-principles calculations 
that describe the nucleation and growth of O-deficient C clusters in Si02 under oxidation conditions, 
fed by the production of CO at the advancing interface, and their gradual dissolution by the supply 
of O under reoxidation conditions. We predict that both CO and CO2 are released during both 
steps. 
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The most significant property of semiconductors is 
their ability to sustain heterogeneous n-type and p-type 
doping. This property, however, is eroded by high tem- 
peratures and high voltages that cause intrinsic excita- 
tion of electron-hole pairs across the band gap. As a re- 
sult, semiconductors with significantly larger band gaps 
than silicon have been investigated as candidates for elec- 
tronic devices suitable for high temperatures and high 
voltages. Silicon carbide is a particularly attractive can- 
didate because its native oxide is Si02 which works so 
well as a dielectric in Si-based microelectronics. The 
presence of a third element, however, namely C, results 
in a wide range of phenomena that do not occur in the 
Si-Si02 system. In particular, oxidation of SiC entails 
the production of CO which effuses through the oxide 
1^,1). Afanas'ev et al. have suggested that carbon clus- 
ters at and near the interface form during oxidation [pLQ , 
but the structure and dynamics of these clusters has not 
been established. Lipkin and Palmour found that, after 
oxidation, a 'reoxidation' step is necessary to produce 
high-quality oxides and SiC-Si02 interfaces P^|. During 
this step, oxygen is supplied as during oxidation, but the 
temperature is lowered so that no further oxidation takes 
place. In contrast, post-oxidation heat treatment with- 
out the supply of O leads to an increase of charged defects 
in the oxide 0. It is believed that the 'reoxidation' step 
cleanses the interface and bulk Si02 of residual carbon 
H . Duscher et al. recently presented direct experimental 
evidence for the existence of carbon in as-grown samples 
and its removal after reoxidation 

The nucleation and growth of impurity clusters in 
semiconductors is a generic problem for which total- 
energy calculations are well suited to provide detailed 
information. In this Letter we present the results of 
extensive first-principles density-functional calculations 
that allow us to give a detailed account of the nucleation 
and growth of O-deficient carbon clusters in Si02 during 
oxidation conditions and their dissolution during reoxi- 
dation conditions. Basically, a CO molecule, generated 



at the advancing interface and diffusing through the ox- 
ide, can bind weakly to an O site in the Si02 network. A 
second CO molecule, however, can bind to the first and 
the new complex is very stable. Additional CO molecules 
can extend the cluster. The process is helped by a passing 
CO that takes an O atom from the cluster and effuses as 
CO2. The net result is O-deficient clusters. During 're- 
oxidation' at reduced temperatures, no further oxidation 
of the SiC substrate occurs so that no CO is produced. 
Instead, the supply of O atoms is used to 'oxidize the 
C clusters'. The process leads to the gradual dissolution 
of the C clusters through the formation of both CO and 
CO2. These results account for the available data and 
lead to a prediction that can be tested, namely the effu- 
sion of CO2 during both oxidation and reoxidation. Note 
in particular that, during oxidation, there is simultaneous 
reoxidation, accounting for the fact that the concentra- 
tion of C clusters is not excessive and that substantial 
amounts of CO are in fact released. In contrast, during 
reoxidation, only dissolution of the clusters occurs. 

The objective of the calculations was to determine one 
or more sequences of steps that lead to the formation of 
O-deficient C clusters by the aggregation of CO molecules 
and simultaneous release of O. Thus, sequences of stable 
configurations were determined. Determination of the 
barriers for the formation of theses configurations is a 
laborious task because of the multitude of atoms that are 
involved. Test calculations were carried out in several 
cases by assuming specific pathways. Typical barriers 
were smaller than 1 eV. Since the calculations represent 
upper bounds and the relevant temperatures are above 
1000°C, the processes identified in this paper do in fact 
occur. 

The calculations were performed within density func- 
tional theory, using the pseudopotential method and a 
plane wave basis set. The exchange-correlation effects 
were treated with local density approximation in the form 
of Ceperley and Alder M . We adopted the Vanderbilt ul- 



trasoft pscudopotentials |10|. A plane wave energy cutoff 
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of 40 Ry was used. All the structures were fully relaxed. 
Similar to earlier work , the calculations for the struc- 
ture relaxation were conducted with one k point at (0.5, 
0.5, 0.5), while the densities of states (DOS) were calcu- 
lated with 14 inequivalent k points. 

A supercell (9.60811 x 10.93795 x 10.62372 A 3 ) contain- 
ing 72 atoms, generated by Tu et al. [Q with a Monte 
Carlo bond-switching method Jl3| , was used to model the 
amorphous SiC>2 . This model corresponds to a density of 
2.15 g/cm 3 for the amorphous SiC>2, close to the experi- 
mental value of 2.20 g/cm 3 at the standard pressure. The 
calculated DOS was found to be in good agreement with 
the experimental photoemission spectra |14[| and previous 
first-principles calculations Jl5|,|l6|], which used different 
structure models to simulate the amorphous Si02- The 
Si02 band gap was calculated to be 5.6 eV, the same as 
in earlier calculations fl6[| . 

CO molecules in a perfect Si02 network. We assume 
that the advancing interface during oxidation emits CO 
molecules, as found in Refs. jjj and (2), and explore the 
possible reactions that CO molecules can undergo in the 
Si02 matrix. A perfectly bonded network without any 
defects is first considered. We find that a single CO 
molecule can bind to a network O atom. The bound con- 
figuration (hereafter [CO]) is shown in Fig. 1(a). There 
is a Si-C-O-Si bridge with the extra O atom attached to 
the C atom so that the integrity of the CO molecule is 
preserved. The bondlength of the CO molecule is now 
1.21 A, 0.07 A longer than the bondlength of a free 
CO molecule at an open interstitial site in Si02 (here- 
after CO( /ree )). While this C-0 bond is more likely a 
double-bond, another C-0 bond between the carbon and 
the original bulk oxygen atom appears to be a single- 
bond (with a bondlength of 1.36 A). The total energy of 
the [CO] configuration is 0.23 eV lower than that of the 
CO( /ree ). This is a rather small binding energy, suggest- 
ing that both structures can exist. The energy cost for 
the release of a CO molecule from Si02 is calculated to 
be 0.18 eV. Hence, CO molecules at the open interstitial 
sites can easily escape from Si02, as found by experi- 
ments jjj. 

The binding energy of a CO molecule in the [CO] 
configuration, though small, is sufficient to result in a 
substantial fraction of CO molecules being bound at 
any given time [about 9 out of 10, as estimated from 
exp(-AE/kT) with AE = 0.23 eV and kT = 0.1 eV]. 
When a diffusing CO molecule encounters a [CO], it can 
be attached to form a two-molecule complex, [C2O2], via 
the reaction 

[CO] + CO (free) -> [C 2 2 ]. (1) 

The energy gain is 1.6 eV, which suggests strong binding. 
The [C2O2] configuration is shown in Fig. 1(b). The 
complex forms a bridge across a network ring with the 
two halves of the bridge being on different planes. 



When a third CO molecule gets near the [C2O2] struc- 
ture, instead of binding, it takes one of the oxygen atoms 
from the cluster and forms either a CO2 molecule at an 
open interstitial site (C0 2 ^ >ee ' ) configuration) or a CO2- 
like bonded configuration ([CO2], shown in Fig. 1(c)). 
The [CO2] configuration has essentially the same total en- 
ergy as CO^ 66 '' (only 0.04 eV lower, which is within the 
uncertainty of the calculations). The remaining bonded 
configuration ([C2O]) contains two carbon and one ad- 
ditional oxygen atoms, as shown in Fig. 1(d). This 
structure is similar to the [C2O2] configuration except 
for having one less oxygen atom. The above process re- 
sponsible for the formation of CO2 is exothermic with a 
energy gain of 0.8 eV. Such a reaction can be expressed 
as: 

[C2O2] + CO (/ree) -> [C 2 0] + CO { 2 free} + 0.8 eV. (2) 

We conclude that oxidation of SiC results in the release 
of both CO and C02- There are no experimental reports 
about CO2 being released. Special effort would be needed 
to detect it because it is probably present in the chamber 
already. 

As more CO molecules arrive, the [C2O] complex un- 
dergoes additional reactions. First, a CO molecule binds 
to the [C2O] structure which results in the formation of 
a [C3O2] configuration, shown in Fig. 2(a). The complex 
contains one C-0 double bond and two C-0 single bonds 
with bondlengths of 1.15 A and 1.31 A, respectively. The 
bondlengths of the two C-C bonds are 1.31 A and 1.37 
A. The energy gain is again very large, 2.0 eV. 

The [C3O2] configuration is a very stable defect. The 
arrival of yet another CO molecule does not lead to fur- 
ther growth. Though another CO can be attached to 
[C3O2], there is an energy cost of 0.54 eV, which makes 
the reaction highly unlikely. Similarly, a passing CO 
molecule cannot take away an O atom from [C3O2] and 
become C02- The process again involves energy cost 
(0.22 eV). If, however, by coincidence, two CO molecules 
happen to arrive at a [C3O2] complex at the same time 
(a rare event), the reaction 

[C 3 2 ] + 2CO ifree) -> [C4O2] + CO [ 2 Sree) (3) 

can occur with an energy gain of 1.3 eV. The resulting 
[C4O2] complex is shown in Fig. 2(b). Once it forms, 
this complex can again lose an O atom to a passing CO, 
becoming a [C4O] complex, shown in Fig. 2(c). The gain 
in energy is 1.4 eV. 

The [C4O] is unstable, however. Upon the arrival of a 
CO molecule, the following reaction occurs 

[C4O] + CO Uree) -» [C3O-2] + 2CO ( 2 free) + 1.1 eV (4) 

with an energy gain of 1.1 eV. Two CO2 molecules are 
released and the residual defect [CsO_2] consists of three 
C atoms in the space of two O vacancies [Fig. 2(d)]. This 
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last defect is stable against the arrival of one or two CO 
molecules. 

In summary, we have shown that through a process 
of complex reactions CO molecules passing through S1O2 
will form oxygen-deficient carbon complexes through the 
emission of CO2 molecules. The most stable defects are 
[C3O2] and [C30_2]- The formation of the latter, how- 
ever, necessitates an intermediate step that requires the 
simultaneous arrival of two CO molecules. 

CO molecules in Si02 with O vacancies. When oxy- 
gen vacancies exist in bulk Si02, a CO molecule fills a 
vacancy by forming two Si-C bonds. The resulting con- 
figuration is shown in Fig. 1(e) (thereafter as [C]). The 
C-0 bond is a double-bond with a bondlcngth of 1.23 A. 
The binding energy of the CO molecule in the [C] struc- 
ture (relative to the molecule in the CO^ ree ^ structure) 
is 2.2 eV. However, this structure is not stable. When an 
additional CO molecule arrives, the [C] structure trans- 
forms into the [C2O] structure [Fig. 1(d)]. The energy 
gain is 3.0 eV. 

Reoxidation. We now examine the reactions that can 
occur during reoxidation when no new CO molecules are 
produced but oxygen is available. 

First, defects containing a single carbon atom such 
as [C] and [CO] can be removed in the form of CO2 
molecules through the following processes: 

[C] + O -» [CO] + 8.8 eV, (5) 
[CO] + 0^ CO ( 2 free) + 8.1 eV. (6) 

Second, clusters of two carbon atoms, ie., those in the 
[C2O] and [C2O2] configurations can be dissociated and 
released in the form of CO and CO2 molecules as follows, 

[C 2 0] + 0^ [C2O2] + 7.6 eV, (7) 
[C2O2] + O -> CO (/ree) + CO ( 2 free) + 6.6 eV. (8) 

Similarly, the three-carbon clusters will be removed 
through the following reactions: 

[C3O-2] + ° -> [C3O-1] + 6.7 eV, (9) 

[C3O-1] + O -> [C 3 ] + 6.8 eV, (10) 

[C 3 ] + 0^ [C 3 0] + 7.5 eV, (11) 

[C 3 0] +0^ [C 3 2 ] + 8.6 eV, (12) 

[C3O2] + O -> [C2O] + CO[ Sree) + 6.3 cV. (13) 

Finally, the clusters with four carbon atoms will lose 
one carbon atom and reduce to three-carbon clusters as 
follows, 

[C4O] + O -» [C4O2] + 7.0 eV, (14) 
[C4O2] + O -> [C 3 2 ] + CO {free) + 7.0 eV. (15) 

The C3 clusters will then be removed with further reox- 
idation, as shown in Equs. (13), (7), and (8). 



In summary, we have shown that the CO molecules 
that are generated during oxidation of SiC diffuse 
through the network and can form a large variety of com- 
plexes containing one, two, three, and four C atoms. One 
key result is that these defects are oxygen deficient, as 
O atoms are removed by combining with CO molecules 
forming CO2 molecules. This phenomenon provides a 
natural explanation of the 'reoxidation' step since we also 
showed that a supply of oxygen can dissolve the C clus- 
ters. 
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FIG. 1. Schematics of several species containing one and 
two C atoms: (a) [CO]; (b) [C 2 2 ]; (c) [C0 2 ]; (d) [C 2 0]; 
and (e) [C] configurations. The dark circles denote C atoms. 
The larger and smaller grey circles are for Si and O atoms, 
respectively. 

FIG. 2. The minimum-energy structures corresponding to 
the configurations containing clusters with three and four C 
atoms in Si0 2 : (a) [C3O2]; (b) [C4O2]; (c) [C4O]; and (d) 
[C 3 0_ 2 ]. 
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